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Abstract 

In this paper we calculate the Bogoliubov coefficients and the energy density 
of the stochastic gravitational wave background for a universe that undergoes 
inflation followed by radiation domination and matter domination, using a 
formalism that gives the Bogoliubov coefficients as continous functions of 
time. By making a reasonable assumption for the equation of state during 
reheating, we obtain in a natural way the expected high frequency cutoff in 
the spectral energy density. 
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1 Introduction 



In the usual framework for the calculation of the stochastic grav- 
itational wave background spectrum |I[ it is assumed that the transition 
between two consecutive epochs in the evolution of the universe is sudden, 
requiring only the continuity of the metric and its first derivative at the tran- 
sition point (when the unperturbed space-time manifold is the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model, this reduces to the continuity 
of the scale factor and its first derivative). Matching the wave function for 
the perturbation and its first derivative at both sides of the transition, we 
obtain the Bogoliubov coefficients relating the initial and the final states. 
This sudden transition approximation gives correct results, except for very 
large frequencies where it gives an overproduction of gravitons that causes 
the energy density to diverge 0. This overproduction can be traced back to 
the fact that in this approximation, the second derivative of the scale factor 
and therefore the scalar curvature is discontinuous ||. To solve this problem, 
a cutoff is introduced for frequencies / larger than the rate of expansion f t 
at the time the transition occurs, by invoking general adiabatic considera- 
tions [fl, p. Loosely speaking, the modes with f ^ ft will not notice that 
the universe is expanding and hence will be exponentially suppressed; this 
suppression is well simulated by the cutoff. But for frequencies near f t the 
problem cannot be correctly solved by the use of a sharp cutoff. 

Contrary to what happens with density perturbations where the high 
frequency part of the spectrum can be easily distorted, the high frequency 
part of the gravitational wave spectrum will remain unchanged during the 
evolution of the universe due to the fact that gravitational waves decouple 
very early. Thus, if our goal is the comparison of our spectra either with 
data from direct observations || or Cosmic Microwave Background Radiation 
measurements (CMBR) ||, more attention should be given to the high 
frequency part of the gravitational wave background spectrum. 

In this paper we calculate the exact Bogoliubov coefficients for the 
production of the gravitational wave background using a method developed 
in ||. This method does not rely on the sudden transition approximation 
and can be used to obtain the Bogoliubov coefficients for the high frequency 
modes without any approximation. Making a reasonable assumption about 
the equation of state during reheating and using a scale factor that interpo- 
lates smoothly between a radiation dominated universe and a matter domi- 
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nated universe we obtain in a natural way the high frequency cutoff. 

The outline of this paper is as follows: in section II we summarize 
our method for the calculation of the Bogoliubov coefficients; in section III 
we calculate the exact Bogoliubov coefficients for a universe that undergoes 
inflation, reheats and evolves through a radiation dominated epoch followed 
by a matter dominated epoch; we end the paper with our conclusions in 
section IV. 



2 The method 

For simplicity we assume that the unperturbed space-time manifold 
is the flat FLRW model. The perturbed metric can be expressed in the 
synchronous gauge as 

ds 2 = a 2 (r) [-dr 2 + (<% + ^(t, x)) dx i dx j ] (1) 

where r is the conformal time, a(r) is the scale factor and hij(r, x) is the 
perturbation (latin indices run from 1 to 3). Expanding hij in plane waves 
we get 



P =i J (2vr)- i / 2 a(T 



d 3 k 

~)V2k 

a p (k, T)£jj(k,p)e lk ' x £(A;, r) + her m. con j. (2) 

where x is the spatial coordinates three-vector, k is the comoving wave- 
number three- vector, k — |k| = = tua, p runs over the two possible 
polarizations of the gravitational waves, Eij(k,p) is the polarization tensor, 
a p (k, r) the annihilation operator and r) the mode function for the grav- 
itational waves obeying the equation 

e + (k 2 - -)£ = (3) 

(X 

where the ' denotes the derivative with respect to the conformal time. The 
annihilation and creation operators can be expressed in terms of the ini- 
tial time annihilation and creation operators A p (k) and Ap(k) through the 
Bogoliubov coefficients a(k,r) and (3(k,r) 

a(k, r) = a{k, r)A(k) + (3*{k, r)A\k) (4) 
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(we have dropped the polarization subscripts for clarity) where the a and (3 
must satisfy 

\a\ 2 -W = l. (5) 

In ref. [f| it has been shown that the Bogoliubov coefficients obey the 
set of coupled ODE's 



a\r) = 4( a(r)+ ^ (T)e2 ^ TO) )l^ (6) 
(3'(r) = -±(p( T )+a(T)e-™^)^ (7) 

where tq is an arbitrary constant. These equations are equivalent to Parker's 
integral equations for a and f3 provided we make for his W(k, r) the ansatz 
W(k, t) = k. Substituting 

a = -{X + Y)e ik{T - T " ) (8) 
(3 = -{X - Y)e~ ik{T - To) (9) 
with X = X(k, r) and Y = Y(k, r) our system takes the form 

X"+(k 2 -^-)X = (10) 

Y = % -X' . (11) 
k 

When the scale factor has the power-law form, the solution of equation (|i0|) 
can be expressed in terms of the Hankel functions |T0[ . Once we have solved 
equation fllCf ) for the entire period under investigation, a and (5 can be easily 
calculated from equations (pTT|) , (||) and (||). To obtain the final Bogoliubov 
coefficients we must now perform one more Bogoliubov transformation to the 
mode functions appropriate to the present state of the universe as explained 
in ref. 0. Denoting the readily calculable Bogoliubov coefficients corre- 
sponding to the final transformation by ct/ t and j3f t , the final Bogoliubov 
coefficients ap and f3p will then be given by 

a.p = aaf t + /3f3f t (12) 
= pa} t + ap ft . (13) 
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With the formalism described above no sudden transition approxima- 
tion is needed and the high frequency part of the graviton spectrum can be 
obtained exactly, provided we can make some reasonable assumption about 
the behaviour of the scale factor during the transitions between two consec- 
utive epochs. 



3 The Bogoliubov coefficients 

We will consider a Universe that undergoes inflation, reheating and 
enters an epoch where dynamics is governed by a mixture of radiation and 
dust. The scale factor has the form 



a(r) = < 



TJT^ V T < T I 

H(t\-t) — 1 

a R (T) Tj < T < T R (14) 

a eq (r 2 + 2r) r > r R 



where Or(t) is an yet undetermined function of r that will smoothly inter- 
polate between inflation and the next epoch, a eq is the scale factor at the 
time of equality between the densities of radiation and dust r eq and H (the 
Hubble constant during inflation) and T\ will be determined later. The so- 
lution is normalized in such a way that r eq = \pl — 1. The scale factor for 
r > tr is a solution of the Friedmann equation for a universe filled with a 
non- interacting mixture of radiation and dust. This assumption provides a 
smooth transition from a radiation dominated universe to a dust dominated 
universe and is the simplest possible improvement over the usual treatment, 
certainly giving a better approximation near r eq ; for r <C r eq the scale factor 
reduces to the one for a radiation dominated universe a oc r and for r r eq 
we recover the scale factor for a dust dominated universe a oc r 2 . 

We still have to do some assumptions about the scale factor during re- 
heating. A reasonable way to procede is to assume an independent equation 
of state for each of the components of the fluid filling the universe and inte- 
grating the Friedmann and energy density conservation equations to obtain 
or(t). A simple but not unique choice for the equation of state is 

Pi (r) = \(r) Pi (r) (i = l,2) (15) 
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where the Aj must satisfy the limiting values 



Ai(r r ) = A 2 (rj) = -1 

Ai(rfi) = | (16) 
A 2 (r R ) = 

as appropriate to the inflationary and radiation+dust stages respectively. 
Of course, our choice of equation of state does not intend to be more than 
illustrative of the method we are using. From the present values of p ra d and 
Pdust we obtain Pi(r R ) 



Pl(rR) = Prad(r R ) = prad{r p ) -p-^r (17) 

P2(r R ) = pdust{r R ) = pdust{r P ) ( 18 ) 

where r p is the present conformal time. 
Denoting by f(r) the function 

/(-) = { < / . T "n (19) 

J v ' \ exp(-ci/T) r > v ' 

we choose \ as 

3f(T-Tj) + f(T R -T) f{ T - Tl ) + f{ TR - T ) 1 1 

A2 " "/(r-r / ) + /(^-r) (21) 

where the constant ci, can be used to control the rate at which the transition 
occurs. It can be easily seen that these A's obey conditions (|i6|). As fig. 1 
shows, these are very smooth functions. 

We can finally write down the equations that determine a R , p\ and p2 



a'\ 2 8ttG 



a 2 J 3c 2 



(P1+P2) (22) 



3-pi(l + Ai) (^ = 1,2) (23) 
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with the X{ given by (|20| ) and These equations can be easily solved 

numerically If we want our model to be consistent with the present day 
observed universe then we must evolve back in time the present condition of 
the universe to obtain the state of the universe at r# and then we integrate 
equations ( p2|) and (|23|) backwards in time up to 77. This procedure will give 
us a scale factor that is compatible with the observed Universe. Once the 
numerical integration is done we can easily find H and 77 appearing in eq. 

a 2 R (n) o fl (rj) 
# = / 7 x , r i = , / x + T i I 24 ) 

In our calculations we used for the redshifts at r eq , tr and 77 the values 
1 + Z eq = 2.38 x 10 4 , 1 + Z R ; = 3 x 10 28 and 1 + Z, = 1.6 x 10 29 . 

Knowing a(r) for the entire period under consideration, we are now 
prepared to calculate the Bogoliubov coefficients. The initial conditions for 

a and (3 are 

H-oo)| = l; |/3(-oc)| = (25) 
which translate for the X and Y variables into 

\X(-oo)\ = \Y(-oo)\ = 1 (26) 

Equation ([10]) can be integrated analyticaly up to r = 77. The solution that 
satisifes (|26|) gives at 77 

= '"'"'fe) (27) 
y(T '» = e - iiT " n ' i 1 + - W^r) (28) 

At this point we must precede with the integration numerically using (p7|) 
and ( P8| ) as initial conditions. After recovering a and /3 from (Rp and (^) 
we must perform the final Bogoliubov transformation ( |I2] ) and (0). To 
do this, we notice first that in spite of our final state being one of mixed 
radiation and dust, as t p ^> r eq the contribution of the radiation component 
is negligible (we have seen at the beginning of section III that for r ^> r eq 
we recover the scale factor for a matter dominated universe a oc r 2 ) and we 
can consider that the final state in our model is one of matter domination. 
As explained fully in ref. |J the final Bogoliubov coefficients must be those 
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appropriate to the present day matter era mode functions and this requires 
a final Bogoliubov transformation to the matter era mode functions with the 
corresponding Bogoliubov coefficients aft and (3ft in equations ( |T2"D and (|13|) 
given by 



«/* = e tkT ° ! + (29) 




-jfe(2r p -To) 

= • (30) 

Figures 2 shows our final results expressed as the energy density 
P{oj) = dE/duj = hu! 3 \P\ 2 /ir 2 c 3 . It is clear from this figure that with the 
approach developed in this paper it is possible to obtain, in a natural way 
the high frequency cutoff. 

For the sake of comparison we also calculate the Bogoliubov coeffi- 
cients in the sudden transition approximation for a universe that undergoes 
inflation, radiation domination and matter domination. Assuming that in 
this case the transition between inflation and radiation occurs at Tj r and 
denoting the Bogoliubov coefficients associated with the transition from in- 
flation to radiation domination and from radiation to matter domination 
respectively with the subscripts r and d we have 

- — f 1 + W, ~ (31) 

* = mr (32) 

= <=""■' ( l + 7^-^-^) (33) 



2kr eg 8(kr ( 



eq ) 



The final Bogoliubov coefficients are then given by fl|].|| 

{cx r k m i n <C k ^ k r (351 

a r a d + (3 r (3* d k r < k < k cut 



and 



n J Pr k m i n <C k ^ k r ( Qfi\ 

P ~ \ (3 r a* d + a r (3 d k r < k < k cut 
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with 



k min = 2ira(r p ) H(t p ) (37) 
k r = 2Tra(r eq ) H(r eq ) (38) 

To fix the value of k cut we notice that, acording to general adiabatic argu- 
ments, the characteristic frequency for the cutoff is given by u> cu t = 2 n At- 1 
where At is the (physical) time scale for the transition, usualy taken to be 
H^Tir)' 1 |1(J], fll!] . However, as these authors point, the exact value of At 
depends on the details of the transition. In our model At is known and we 
have At w a(r ir )Ar which gives 

hut = ^- • (39) 

At this point we still have to determine n r . Because in our model 
we have introduced one intermediate epoch between inflation and radiation 
there is, when comparing the two models, an ambiguity in the choice of Tj r . 
To overcome this ambiguity we choose r ir in such a way that the number of 
gravitons at the beginning of the radiation dominated era is the same in the 
two models. 

Figure 3 shows the comparison of the results obtained with our method 
and the one using the sudden transition approximation. We can see in Fig. 
3a that the method discussed in this paper predicts a larger P for the low 
frequency gravitons than is the case with the usual method of calculation. 
This can be easily explained by the fact that during the radiation dominated 
epoch, and if we assume only radiation to be present, there is no graviton 
production due to the well known conformal invariance of the equations 
PI , [[0|. In our model the invariance is broken because we also have dust 
mixed with the radiation and we thus have an additional contribution to P. 
No matter how small the contribution of the dust component becomes, the 
conformal invariance will always be broken, as can be seen by the fact that 
a" = a eq in conformal time. We should also notice (fig. 3b) that the point 
where the high frequency suppression starts to take place, depends on the 
time it takes for the transition between inflation and radiation domination, 
as expected from adiabatic arguments [|J ||. The dots in Fig. 3b represent 
the energy density obtained with the sudden transition approximation and 
k cu t given by ( |3"§D with At = 6. x 10~ 26 , corresponding to the width of the 
transition in our model with C\ = 11.. Comparing the two curves, we can 
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see that, for frequencies near the suppression, the energy density obtained 
by our exact calculation is one order of magnitude bellow the result obtained 
with the sudden transition approximation. 



4 Conclusions 

We have shown in this paper how the high frequency cutoff in the 
energy density of the stochastic gravitational wave background could be ob- 
tained in a natural way. Although our results are not surprising, the calcu- 
lation done in this paper is exact, even in the high frequency region of the 
spectrum. We showed that for frequncies slightly below the high frequency 
cutoff, the sudden transition approximation predicts an overproduction of 
gravitational waves. Differences to previous calculations were also obtained 
in the region of frequencies corresponding to the transition between the ra- 
diation dominated and the matter dominated epochs. 



After this work was completed we came across a preprint by Koranda and 



Allen [14] where the epoch following inflation is also treated as a mixture of 



radiation and dust. 
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FIGURE CAPTIONS 



Figure 1 The A's given by equations ( p0[) (a)and ( pTCf ) (b) as functions of the 
conformal time r during the transition from inflation to the radiation+dust 
era. We can see that the width of the transition can be controlled with the 
parameter c\. Larger values of c\ correspond to faster transitions. We used 
Tj = 3.9585 x 1(T 27 and r R = 3.9585 x 1(T 25 in all the figures. 

Figure 2 The energy density of the gravitational waves produced during 
the expansion of the universe obtained with the method developed in this 
paper. The steeper part of the curve (/ < 10 _16 s _1 ) corresponds to gravitons 
produced during the radiation+dust epoch and the region with constant slope 
is due to graviton production during inflation. The sharp suppression is due 
to the finite width of the transition from inflation to radiation domination. 
This figure was obtained with c\ = 11. 

Figure 3 For low frequencies, (/ < 10 _16 s _1 ) the contribution for P of the 
gravitons produced during the radiation+dust epoch is larger in our model 
than in the sudden transition approximation. Due to the presence of a small 
contribution of dust in the radiation dominated era, the conformal invariance 
of the equations is broken giving rise to more graviton production (a). As 
it was expected, the point where the high frequency suppression appears 
depends on the time it takes for the transition from inflation to radiation 
domination. In our model the width of the transition can be varied with 
ci (b). The dots represent the results obtained with the sudden transition 
approximation; we introduced a high frequency cutoff k cut = 2ti/At with 
At = 6. x 10 -26 corresponding to the width of the transition for the case 
ci = 11. (See Fig. 1). 
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